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Abstract 

If the scale dependence of a Yukawa matrix is assumed to be de- 
termined entirely by the dominant 33-element, then the renormalization 
group equations can be expressed in terms of two separate equations: a 
differential equation for the 33-coupling, and, an algebraic equation for 
the scale-independent 3x3 matrix that is found to have only two non- 
trivial, hierarchical, solutions with eigenvalues (0,0,1) and (0,1,1). The 
mass matrices are constructed from these solutions by rotating them first 
by the experimentally known mixing matrices-the CKM for quarks and 
charged leptons, and the CKM-analog for the seesaw generated Majorana 
neutrinos-and then incorporating the appropriate texture zeros. A uni- 
form, hierarchical, description for the mass matrices of quarks and leptons 
is thus achieved, in terms of the mixing paramters, that give mass eigen- 
values consistent with experiments as well as reproduce the input mixing 
angles. Inverted hierarchy in neutrinos is also discussed. Only a single 
scale 10" GeV)ioT the seesaw neutrinos is involved rather than their 
mass distribution. No new particles are otherwise invoked. 



1 



I. Introduction 

The Yukawa matrices which describe the mass values in the quark sector of 
the standard model show a pronounced hierarchical pattern, as do the mixing 
angles of the CKM matrix [1, 2], 

In the lepton sector the charged lepton masses show a similar pattern of 
hierarchy. The neutrinos, on the other hand, are massless in the minimal ver- 
sion of the standard model, as no right-handed neutrinos arc assumed to exist. 
The neutrino oscillation data [3,4,5,6,7,8] indicate, however, that neutrinos have 
masses which are extremely small and have a pattern which can be consistent 
with hierarchy or, alternatively, inverted hierarchy [9]. 

The mixing angles, the CKM-analog [10], for leptons for a diagonal charged 
lepton mass matrix are quite different. While for CKM, all three rotation angles 
are small, here two of the angles arc foimd to be quite large[3,4,5,6,7,8]. 

We would like to address all these facts together and find a uniform, underly- 
ing explanation by solving the renormalization group equations in MSSM (mini- 
mal supcrsymmctric standard model) [11,12] supplemented by two assumptions 
which reflect physics beyond the standard model. We assume the well known 
seesaw mechanism [13] for generating the small (Major ana) neutrino masses 
in which the standard model left-handed neutrinos couple to large mass right- 
handed neutrinos. And we take account of the fact that texture zeros may 
be present in the Yukawa matrices of quarks, charged leptons and neutrinos 
[1,11,14,15,16,17,18,19]. 

We begin first in, sections II and III, with discussing the renormalization 
group equations (RGE) of the Yukawa matrices. The RGEs, to one loop in 
MSSM, for up-quark (f/) and down-quark (D) Yukawa matrices are given as 
follows [11] 



dU _ _J_ 
H ~ 167r2 



agf + 3UU+ + DD+ + Tr{3UU+)l 



U (1.1) 



dD 

Itt 



167r2 



- ^ c^flf + 3DD+ + UU+ + Tr{3DD+)l 



D 



(1.2) 



where t = ln{iJL /IGev), ji being the energy variable, Cj and c[ are known con- 
stants and g[s are the gauge couplings. 

In the lepton sector, we consider the charged lepton Yukawa matrix, i?, and, 
for the neutrinos we assume a seesaw mechanism, as mentioned earlier, and 
consider the (Dirac) neutrino matrix, N , that couples the standard model left- 
handed neutrinos to the large-mass right-handed seesaw neutrinos given by the 
mass distribution, M_r.To one loup in MSSM the RGEs are given by [11,12] 



Hit 



167r2 



~Y^digf + 3EE+ + NN+ + Tr{EE+)l+Tr{3DD+)l 



E 



(1.3) 
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dN 
lit 



167r2 



■ ^ d'igf + 3NN+ + EE+ + Tr{NN+)l+Tr{3UU+)l 



U 



(1.4) 



where djand d' arc known constants. The Majorana neutrino mass matrix, k, is 
then given by the seesaw formula [12,13] 



(1.5) 



Renormalization group equations, used as a tool to determine the properties 
of the quark Yukawa matrices have been considered previously, specfically for 
the case in which the top coupling is assumed dominanat[1.14,20]. Recently, the 
RGEs have also been used to reconcile with the neutrino data and to investigate 
the differences between CKM and the CKM-analog mixing angles [10,12,21] 

Let us consider the U matix given by (1.1). Assuming the 33-matrix element 
to be dominant and equal to the top-quark Yukawa coupling, At, we obtain the 
following equation ignoring D and the g'^s as they are small compared to At 

dXt 3At , , 



The solution is given by, 

(1.7) 

where we have taken At {to) = Xot at a convenient scale parameter to- 

The question we now wish to explore is this: just as the magnitude of Aj dom- 
inates the U-matrix, are there solutions to (1.1) such that the scale-dependence. 
At {t) given by (1.7) also describes the scale-dependence of the entire U- matrix 
? In other words, does equation (1.1) allow a factorizable solution of the type 

U = UoXt (i) (1.8) 



At {t) = A, 



Of 



3 . 



t) A 



that the t-dependence, as given by Af (t) , can be factored out, leaving behind a 
matrix of coefficients Uq that is independent of the scale, t. Similarly, we wish 
to explore the possibility that the solutions to (1.2), (1.3) and (1.4) can also be 
expressed in a factorizable form with the scale-dependence residing entirely in 
the dominant matrix element of each, given by the b-quark, r— lepton, and the 
largest mass (Dirac) neutrino, ms (or mifor inverted hierarchy) e.g. 

D = DoXb (t) (1.9a) 
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E = EoXr {t) 



(1.9b) 



N = NoXav (i) 



(1.9c) 



Such solutions, indeed, exist as we will discuss below, in which the RGEs 
split into two equations, one a differential equation for the dominant (scale- 
dependent) Yukawa coupling and the other an algebraic equation for the (scale- 
independent) matrix of the coefRcients. It is found that the scale-independent 
matrices (e.g. Uo,Dq etc.) can be classified in terms of their eigenvalues in 
simple, diagonal -"primordial"- forms. 

In deriving these results we assume U. D. N and k to be real and symmetric. 

There are, actually only two non-trivial solutions for the scale-independent 
matrices : one we designate as "hierarchical", and the other "semi- hierarchical", 
of the forms 
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(1.10) 



respectively. In particular, Uq is identified as the former, Do as the mixture of 
the two and similarly for the leptons. A general, real, symmetic, solution can be 
obtained from this "primordial" matrix by rotating it through arbitrary angles 
by a unitary matrix. That is, if Udiag and Ddiag are related to the "primordial" 



matrices then one can write the general matrices as 

Uo = VUdiagV^ 



(1.11) 



Do = VDdianV^ 



(1.12) 



where y is a unitary matrix. A similar situation exits for the leptons. 

Thus the Yukawa matrix elements will depend directly on the rotation pa- 
rameters of V. 

In section IV, confining first to the quark sector, what we find most re- 
markable is that, if we assume the unitary matrix V to be the same as the 
CKM matrix, we obtain Yukawa matrices Uq and Dq that are exactly what one 
generally expects when expressed in powers of the Cabibbo parameter, A [1]. 
Namely, 



A'^ A^ 
A^ A^ 



4n 



A 
A^ 
1 



(1.13) 



Similarly for Dq where the hierarchy is found to be less pronounced, just as ex- 
pected [1]. 
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In the Icpton sector, discussed in V, the sohitions for the charged leptons 
exhibit a similar hierarchy through the small angles of the CKM matrix. For the 
neutrinos, we argue that for the mass matrix, k, defined in (1.5) the candidate 
for the appropriate unitary matrix V is the large-angle lepton mixing matrix, 
the CKM-analog [10]. 

We thus show that for each of the four systems of particles, the mass matrix 
can be described in terms of the appropriate mixing parameters through an 
expression that is very simple and transparent. 

Although our model reproduces the expected hierarchical pattern of the 
mass matrix elements, it is, at this stage, incomplete, because, for example, 
the eigenvalues of Uq and Dq in equation (1.11) and (1.12) are still given by 
their "primordial" representatios i.e. (0,0,1) and (0,1,1). Furthermore, the 
mixing matrix turns out to be a unit matrix, since the same matrix, V, that 
diagonalizes Uo, also diagonalizes Do- 

All that is remedied, however, once we incorporate texture zeros, as we do 
in section V, whereby, in Uq and Do, constructed through equations (1.11) and 
(1.12), some of the entries are replaced by zeros. 

Confining ourselves, briefly, to the quark sector we note that the (symmetric) 
Yukawa matrices Uq and Dq that are consistent with experiments and allow the 
maximum number of texture zeros, are found to be of the following type, 



X 
0X0 
X X 



Do 



0X0 
XXX 
X X 



(1.14) 



a total of five texture zeros, where X denotes non-zero entries[16.17]. We assume 
this to be the pattern of the zeros in our case. 

Two things happen as soon as texture zeros are incorporated. First, the 
eigenvalues values of (the newly textured) Uq and Dq will no longer be given 
simply by the primordial values, but will be proportional to the matrix elements 
of V that were involved in (1.11) and (1.12). Second, the mixing matrix will 
not necessarily be a unit matrix since there will now be a mismatch between 
matrices that diagonalize Uo and Dq, as is evident from the structures in (1.14) 

We need then to define two types of mixing matrices '-Vckm ^ckm- 

The "input" mixing matrix, V^!'^\j, represents V which generates Yukawa 
matrices Uq and Do through equations (1.11) and (1.12). The entries in these 
matrices are the experimentally known CKM parameters. 

After texture zeros arc incorporated, according to the structures (1.14), an 
"output" CKM, V^°^*^ is defined which is a mixing matrix between (the newly 
textured) Uq and Do- That is, if Vu and Vd diagonalize the two (newly textured) 
matrices Uq and Dq, respectively, then, 

V^k2, = V^V, (1.15) 
If V^°^2i turns out to be the same as the experimentally known CKM matrix 
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(i.e. if V^'j^ipj = V^*^^) then it is a success of the model, justifying the choice 
of Vq^j^ as the unitary matrix V in (1.11) and (1.12). That is, indeed, found 
to be the case when we carry out our calculations. 

In other words, what was put in to generate the mass matrices, came out, 
in a self-consistent manner, when we calculated the mixing angles. 

At the same time, we also find that the eigenvalues, or rather, the ratios 
of the eigenvalues to the 33-elements are consistent with the experimentally 
determined mass ratios. 

For the leptons, discussed in VI and VII, we take the charged lepton structure 
for E(j to be the same as Dq. And for the neutrinos wo refer to an extensive 
analysis of the texture in the neutrino mass matrix which points to two possible 
types of structures that are consistent with experiments[18]. The so called A- 
type structure given by 
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X 
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or kq 






X 





X 


X 


X 





X 


X 



(1.16) 



gives hierarchical neutrino mass values, and the C-type 
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X 
X 
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X 



X 

X 




(1.17) 



gives an inverted hierarchy [18,19]. 

As in the case of quarks, as we discuss in section V, we have Vq^^ that 
generates charged leptons Yukawa matrices, iio, while for kq) that role, we argue, 
will be played by V^^^g which is the large angle mixing matrix obtained from the 
neutrino oscillation data [3,4,5,6,7,8,10] in the basis where the charged leptons 
are diagonal. And, Vi°^gl will be an appropriate product of the matrix that 
diagonalizes ko, and the one which diagonalizes Eq since the Eq matrix we 
obtain from our prescription will not necessarily be diagonal. 

In other words, the process is the same as for Uq and Dq. We first create Eq 
and Kq by the appropriately designated V. Then we replace some of the entries 
by zeros as given above and examine the newly textured matrices. 

Once again we find that, with the A-type structure, Vi°^gl is the same as 

^large ' sxcept for a Small discrepancy in the the angle S2 for the 1-3 sub matrix. 
The mass eigenvalues are found to be consistent with experiments. 

We discuss the C-type structure for kq separately, in section VIII, since 
it produces an inverted hierarchy which can not be accomodated within the 
framwork of the "primordial" matrices we have pursued which are basically 
hierarchical in nature. Instead we point out that, for this case, the appropriate 
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basis is provided by 



1 





(1.18) 



We find that, as a consequence, V^o°ge essentially identical to Vj^™ g with a 
maximal atmospheric angle mixing (2-3 submatrix), and S2 = 0. The neutrino 
mass eigenvalues are also consistent with experiments. 

II.The RG Equations 

First let us define the following more convenient variable 

^=i^(*o-0 (2.1) 



then, (1.6) and (1.7) can be expressed as 



C^At ^ 1,3 

dx 2 



= o^t (2-2) 



Xtix)=Xot[l + x\l] ^ (2.3) 

Secondly, ignoring the gauge terms, which are small, we re-arrange the re- 
maining terms in (1.1) and (1.2), as follows 

= 1 [UU+ + Tr{UU+)i\ U + ^DD+U (2.4) 



-— = -\DD+ + Tr(DD+)l]D+—UU+D (2.5) 

where the first term, in the square brackets in each, involves "uncoupled" terms, 
and the second involves coupling between U and D. 

In a similar fashion we write, E and N in the following form 



'ir = ~o [^^^^ + Tr{EE+)i\ E + ^NN+E + ^Tr{DD+)lE (2.6) 



= i- [2,NN+ + Tr{NN+)l] N + -^EE+N + ^Tr{UU+)lN (2.7) 
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where, besides the uncoupled and coupling terms, we have a term, in each, which 
contains a trace in the quark sector that basically normalizes (re-scales) the solu- 
tions. 



III. Solutions to the RGEs and the eigenvalues in the quark sector 
The equations involving only the uncoupled terms in both (2.4) and (2.5) 
look like 



dM 1 



[MM+ + Tr{MM+)i\ M 



For a solution of the type (1.8) and (1.9a) 

M = MoXm {x) 



(3.1) 



(3.2) 



Xm {x) representing the dominant 33-matrix element, we have two equations 

(3.3) 



"^^rn jx) _ 1 3 . ^ 



and, 



Mo 



MoM^ + trMoM^ol 



Mo 



(3.4) 



Let us now consider equation (3.4) and assume to be real. We diagonalize 
the equation through unitary matrices Vi and V2 to obtain 



Mdiag = VIM0V2 



(3.5) 



Mdiag = \ [Ml,g+trMl^gl] Mdrag 

we express the above relation in terms of the eigenvalues 

Mdiag 



Ai 

A2 
A3 



then we obtain the following relations 



1 



Ai = -[2A? + Ai (A^ + A^) 



A2=2 N + A2 {Xl + Xl)] 



(3.6) 



(3.7) 

(3.8a) 
(3.8b) 
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A3 = ^[2A^ + A3 (Xl + Xl)] 



(3.8c) 



We assume the eigenvalues to be real and order them in the sequence Ai < 
A2 < A3, and consider only positive values for these fermions. There are then 
the following four possibilities 

(i) All the AiS are zero. This is a trivial solution 

(ii) All the A^s are non-zero in which case they must all be equal. This is 
also a trivial solution if Vi = V2, as is the case for Uo and Dq defined in (1.11) 
and (1.12) and, when Mq becomes a unit matrix 

(iii) Ai = 0, in which case A2 = A3 

(iv) Ai = A2 = 0, in which case A3 = 1 

Therefore, there are only two non- trivial solution, which we describe as fol- 
lows 



(i) "Hierarchical" solution with Ai = A2 = 0, A3 = 1 



M. 



(1) 



diag 



000 
000 
1 



(3.9) 



(3.10) 



(ii) " Semi- hierarchical" solution with Ai = 0, A2 = A3 




1 
1 



(3.11) 



(3.12) 



A. Hierarchical solutions 

In Appendix I a general form of Mq'^ is obtained by using equation (3.4) 
directly. 

Below we obtain Mg^^assumed real and symmetric, from relation (3.5) by 
taking 



Vi=V2 = V 



(3.13) 
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where we will take V as the product of three rotations that diagonalizes, in 
succession, the 2-3, 1-3 and 1-2 sub-matrices by angles 0i,02,and respectively 
i.e. 



C2C3 
-C1S3 - S1S2C3 
S1S3 - C1S2C3 



C2S3 
C1C3 - S1S2S3 
-S1C3 - C1S2S3 



S2 
S1C2 
C1C2 



(3.14) 



where Si and Cj are the corresponding sine and cosine- values. With (3.13) and 
(3.14) we then have the following simple form 



§2 C2S1S2 C1C2S2 

C2.S1.S2 C^sl C1C2S1 

C1C2S2 cic^si cjcl 



(3.15) 



There is thus a direct connection between the mass matrix and the rotation 
parameters. This shows a classsic hierarchy pattern, particularly for the case of 
small angles, in which case one can write (cj ~ 1) 



4 

S1S2 

S2 



S1S2 
Si 



S2 
Si 
1 



(3.16) 



B. Semi-hier£trchical solutions 

Following the same procedure as above, the semi-hierarchical patterns are 
created by writing 



With (3.13) and (3.14), M^'^ in (3.12) becomes, for small angles (cj w 1) 



^3 

S3 

S2 - S3S1 



S3 



-S2S3 



S2 - S3S1 

-S2S3 



(3.17) 



(3.18) 



IV. U and D matrices, and the CKM angles 
A. U-matrix 

If we assume Ab << A* then the contribution of the second term in (2.4) 
given by 

—DD+U (4.1) 
12 ^ ' 

will be of the order j^A^A^ which is much smaller than the first term in (2.4) 

which is of the order^A^ . It can, therefore, be safely neglected. The U-matrix 

will satisfy, to an excellent approximation, the following 

~S " \ [^^^ ^ Tr{UU+)\] U (4.2) 
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which is the same equation as (3.4) for Mq and will have the solutions already 
discussed. 

We now assume that Uq is of the "hierarchical form" so that 



■lag 





1 



and 



Uo = V 





1 



Ft 



If we take V to be the same as the CKM matrix 



(4.3) 



(4.4) 



V = VcKM 

then, as discussed in III A, equations (3.15), wc have 



C2S1S2 C1C2S2 



C2S1S2 



„2„2 
C2S1 

C1C2S2 CiC^Si 



cic|si 

„2 2 



(4.5) 



(4.6) 



where the angles now are the same as given by the CKM matrix , we call it "in- 



(4.7) 



put" CKM, giving us the expression (with 


a R 


^1) 
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S1S2 


S2 




Uo = 


S1S2 


sl 


Sl 






_ S2 


S\ 


1 




with 










T/(»") _ 

^CKM — 


1 




S3 


S2 


-■S3 
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Sl 




_SlS3 - 


S2 


-Sl 


1 



(4.8) 



Let us elaborate further on the above result by expressing the CKM matrix 
in the Wolfenstein representation [2], keeping only the leading terms 



VcKM 



A 

AX^ (1 - p - iv) 



A 
1 

-AX" 

r{in) 



AX-^ {p - irj) 
AX' 
1 



(4.9) 



Comparing this with the expression for Vqj^'j^^ above we can extract si, S2, 
and S3 as follows considering only the magnitudes of the Wolfenstein parameters 
(ignoring phases). 

Sl w AX' (4.10a) 



S2 « aaVp' + »?^ 



(4.10b) 
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S3 ~ A 



(4.10c) 



Substituting the above values in the expression for Uq given by we obtain 



Experimentally [22], one finds, in order of magnitude terms. 



(4.11) 



Awl 



(4.12a) 



therefore, the above expression simpHfies to 



(4.12b) 



A« A^ A^ 
A« A^ A^ 
A* A2 1 



(4.13) 



which is in excellent agreement with the behavior one expects for Uq [1] . 

To recapulate then, after including the x-dependence from (2.3), we have the 
complete expression for U given by 



U 



SlS2 
S2 



SlS2 
s2 

Si 



S2 
Si 
1 



Aot [1 + xAj 



(4.14) 



where Si are the "input" CKM parameters. 

The eigenvalues of Uo are, of course, the same as M^J]^ i.e (0, 0, 1). To obtain 
the correct values, as we stated in the Introduction, one must go beyond the 
standard model which we propose to do in Section V by incorporating texture 
zeros. 

B. D-matrix 

We assume here that D is different from U by the fact that the solution of 
the "uncoupled" part 



-^ = l[DD^+TriDD^)l]D 



(4.15) 



IS now given by M^^\the " semi-hierarchical" solution (3.11) and (3.12) with 
Vi = V2 = V 




1 
1 



12 



However, this can not be a complete solution for the -D-equation given by (2.5) 

= 1 \DD+ + Tr(DD+)l] D + —UU+D 
dx A ^ ■' 12 



since the coupling term neglected in (4.15) 



1 



-UU-D 



(4.16) 



can be very large because At >> Af,. This is in stark contrast to U where the cou- 
pling term was negligible. 

If we consider the x-independent (scale-independent) coupling term in (4.16), 
then we have 



(4.17) 



- - 
12 V 3 



Therefore, 





"0 





0" 








"0 





0" 








"0 





0" 




V 















V 















V 





1 















1 














1 














1 





(4.18) 



12 " ° " 12 V 3 





"0 





0' 




V 





















1 





1 

12 



(4.19) 



So there is a contribution from Mq^-* in the D-matrix that needs to be added 

(2) 

to Mq . One can, to a good approximation, write the differential equation for 
D in (2.5) as 



dx 







"0 





0" 




n 





1 















1 



+ T2V3^*^^ 



"0 





0" 





















1 





= -AgMf + -y^At^A(,M^ 



(1) 



(4.20) 
(4.21) 



We can solve for the 33-elements on both sides of the equation (4.21), as- 
sumed dominated by A^ and At, and obtain 



dXb 
dx 



(4.22) 
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An exact analytical solution of (4.22) is given in Appendix II 
It is important to note that, as far as the 22-elements of the two diagonal 
matrices in (4.20) arc concerned, the second matrix is a "small perturbation" on 
the first because it is, in fact, negligible (zero) compared to the first . In contrast, 
with respect to the 33-matrix elements, it is the first matrix which is a "small 
perturbation" since A5 << A(. Furthermore, we note that the x-dcpcndcnce of 
the 22-diagonal element will not be the same as that of the 33-element, which 
is inconsistent with our assumption that the entire D matix is represented by a 
single function of x. 

However, since the first matrix is extremely small we can write an approxi- 
mate solution for D as follows 



^diag — 




e 
1 



(4.23) 



Do = V 




e 
1 



(4.24) 



e is assumed to be a constant, independent of x, and where in the 33-element 
of Ddiag we have ignored the contribution of e compared to 1. 

To determine e one can (numerically) integrate the first term in (4.22), ^Xf, 
from the knowledge of the analytic expression for Ab given in Appendix II, and 
compare it to the integral of the second term, -^A^ A^, with Aj given in (2.3). 
The ratio of the two will estimate the value of e. If we take the ratio of the 
Yukawa couplings A06 and Aot to be the same as the ratio of the masses then 
we obtain e w .01 which is found to be essentially independent of x. 

Another, rough order of magnitude, estimat(^ of e can also be obtained on 
the basis of the ratio of first to the second term in (2.5) above 



A 



mt 



.025 



(4.25) 



The square root is taken because, effectively, the first term in the equation 

(ID 

(2.5) for — — is proprtional to while the second is proprtional to U^D, so 
dx 

the ratio is w D^. Therefore, the contribution to D will involve the square 
root. The numerical value of e from (4.25) is obtained by taking the ratio of 
the Yukawa couplings to be the same as the ratio of the corresponding masses. 



In our calculations to follow we will take e as an arbitrary parameter keeping 
in mind, however, that its order of magnitude is w 10~^ . 
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To obtain an approximate solution for Ah,wc can ignore the first term in 
(4.22), which is very small, to write, after substituting expression (2.3) for Aj , 



T^A?Ab =-AgJl + a;A^,] A (4.26) 

the solution to which is 



A(, (x) = Ao6 [1 + xA^,] (4.27) 



Ao6 = Afe (to) = Afe (x = 0) 

From (4.24) we can write down the expression for D 



D 



eS3 S2 - esiS3 



es3 
S2 - esiss 



Sl 

1 



Ao6 [l + xAo(] 



(4.28) 



where Sj's are the "input" CKM parameters. We can express s, in terms of A, 
as we did previously when we considered Uq, to obtain 



Da 



eA 



eA A^ 
A^ + e 



eA^ 



A^ - eA^ A^ 



1 



(4.29) 



The estimate for e given by (4.25) is, numerically « A , and, correspondingly 
we have 



Do 



A* A^ A^ 

A^ A^ A^ 
A^ A^ 1 



(4.30) 



This expression shows that the hierarchy in Dq is not as pronounced as it 
was found in U, which is consistent with generally accepted form for Dq in [1] 



Dq 



X' X'- 
x^ x^ x^ 

X^ X^ 1 



(4.31) 



Our results arc in excellent agreement with the above if we ignore the slight 
discrepancy in the 13- and 31- components which, in fact, are quite inconse- 
quential as we will see below when we discuss texture zeros. 

The eigenvalues of Do are given by the eigenvalues of Ddiag, i-c. (0, e, 1). As 
pointed out in the Introduction, to obtain the correct values one must go beyond 
the standard model which we propose to do in Section V by incorporating 
texture zeros. 



15 



V. Mass eigenvalues and CKM matrix in the quark sector with 
texture zeros 

As stated in the Introduction the most acceptable structures with texture 
zeros are the following 





' Q X 








X 







0X0 




, -Do 




X X 


X 




X X 








X 


X 


In terms of our results in section IV we then have 








"0 S2 






" 


eS3 





Uo = 


sj 




£•0 = 


eS3 sl + e 


Si 




S2 1 











1 



(5.1) 



A. CKM matrix 

The matrix elements above involve the "input" CKM parameters Sj already 
defined (4.8) 



(in) 



CKM 



1 S3 S2 

-S3 1 Si 

S1S3 - S2 -Si 1 



(5.2) 



The numerical values of the parameters are known from experiments [22] and 
are given by 



si = (.038 - .044) , S2 = (.0025 - .0048) , S3 = (.219 - .226 ) (5.3) 



The "output" CKM is given by the standard definition 

V^T,!, = V:V, (5.4) 



Vd is the unitary matrix which diagonalizes -Do; and Vu is the unitary matrix 
which diagonalizes Uq. In terms of the rotation angles it is given by (assuming 



T/(o«t) _ 



Si S' 



1^3 



4 



S2 

S'l 

1 



(5.5) 



primes are used to distinguish from the "input" case. 

It is now easy to observe from the expressions in (5.1) that the rotation for 
diagonahzing the 23-sub matrix involves only Dq, and for the 13-sub matrix it 
involves only Uq. Therefore, by the nature of the location of the texture zeros 
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a (complete) mismatch is achieved between the rotation angles for Uq and Do 
for these submatrices. We then have 

s'l = Si, 82 = S2 (5.6) 



mismatch continues for the 12-submatrix, where the rotation only involves Dq 
and we have 

as we discussed in section IV the magnitude of e was estimated to be ~ 10~^ 
which is much larger than sf (w .0016) .Therefore, 

4 « S3 (5.8) 



we conclude that 

^CKM — ^CKM y^^-^J 



B. Mass eigenvlues 

Plugging in the following "input" values 



si = .04, S2 = .004, S3 = .22 (5.10) 



which are within the range of the experimental values given by (5.3) and taking 
e = .05, we obtain the following mass ratios, keeping in mind that the 33-matrix 
elements A* and Af, are normalized to their observed mass values ("exp" means 
experimental values). 

^ = 4.8 X 10-^ (exp « 3 x 10"^) 
mt 

— = 1. 9 X 10"^ (exp w 8 X 10"^) 
mt 



— = 2.4 X 10"^ (exp w 2 X 10"^) 

mb 

— = 5. 5 X 10-2 (exp « 3 X 10"^) (5.11) 
rub 

The agreement with experiments is very good considering the fact that no 
attempt was made to do a detailed numerical analysis to fit all the data. 
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VI. Solutions to the RGEs in the lepton sector 

The RG equations we wish to solve for N and E are given by (2.6) and (2.7). 
We notice that the last term in each equation involves a trace which multiplies 
equally all the matrix elements of N and E. Therefore, it will effectively change 
the scale, though we realize that this is not completely correct since there are 
also non-linear terms present in the equation. We will ignore the last terms in 
(2.6) and (2.7), nevertheless, as a simplifying assumption. We then have the 
following 



^ = — [3EE+ + Tr{EE+)l\ E + —NN+E 



(6.1) 



dN_ 
dx 



[3NN+ + Tr{NN+)i\ N + -^EE+N (6.2) 



solving these equations we want to follow, as closely as possible, the analogy to 
U and D, in particular to identify N with U, and E with D. 

Considering first, as we did for U and D, the uncoupled terms in the square 
brackets. We find that they satisfy an equation very similar to (3.4), except for 
the numerical factors. 



Following the same procedure as before, we write 

L = LqXi {x) 



(6.3) 



(6.4) 



A; {x) represents the dominant 33-matrix element, and, 

1 



3Lo4 + *^^o4l 



Lo 



(6.5) 



(6.6) 



■lag 



(6.7) 



We express the above relation in terms of the eigenvalues 



-I'd; 



lag 



Ai 











A2 











A3 



(6.6 
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then we obtain the foUowing relations 

1 



Ai = -[4A? + Ai (A^ + A^) 



^3 = T [4A^ + A3 (A? + a: 



(6.9a) 
(6.9b) 
(6.9c) 



Once again we have only two non-trivial solutions. 



(i) "Hierarchical" solution with Ai = A2 = 0, A3 = 1 



r(l) 
^diag 



000 
000 
001 



(6.10) 



(6.11) 



(ii) " Semi- hierarchical" solution with Ai = 0, A2 = A3 



r(2) ^ 

diag 



000 
1 
1 



(6.12) 



(6.13) 



We now make the following important identifications: 
(i) As we mentioned earlier, we identify N with U so that we take it to be 
of the hierchical type i.e. 



AT.. - - 






1 



(6.14) 



(ii) We, therefore, take the 33-element of N much larger than that of E, 
which implies that 



A3?) >> A^ 



(6.15) 
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(iii) Wc observe that unlike U, D, E, which arc entirely standard model con- 
stituents, N couples a standard model object, the left-handed neutrino, i^l, to 
a particle not in the standard model framework. This leads us to identify Vi as 
the small angle rotation matrix of the standard model i.e. VckmA^ wc did for 
U, D and E and we assume that V2 represents the large-angle mixing matrix 
of the neutrino data i.e. the analog of the CKM matrix [10] , we have called it 
Viarge- Therefore, 

Vi = VcKM (6.16) 



V2 = Vla 



(6.17) 



from (6.11) 



-(1) 



CKMJ^diag^large " ^CKM 



Vn 





1 



V, 



large 



(6.18) 



N = NoX3v{x) (6.19a) 

Wc can write down the expression for X^^ in analogy to U taking account of 
the differences in the numerical factors in the equation for U in (2.4) and for N 
in (6.2), 

X3v{x)=Xo3,[l+xXls^]~~' (6.19b) 

where, Xo3v = A3„(0). 



(iv) We identify E with D, and write 



e 




(6.20) 



e',as in the case of D, is the mixing parameter between the hierarchical and semi- 
hierarchical matrices. And, of course, since it is a standard model particle, we 
take, as we did for Uq and Dq, 



Vi =¥2 = VcKM 



(6.21) 



Therfore, 



£'0 = VcKM 




e' 
1 



CKM 



(6.22a) 
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E = EoXr{x) 



(6.22b) 



Finally, for the Majorana neutrino mass matrix k, described by the seesaw 
mechanism, as given by (1.5) we have 

k = N'^ [M^^] N (6.23) 

as mentioned in the Introduction, [^'^/^^j is the distribution of the reciprocal of 
the seesaw neutrino masses. We then have from (1.5) and (6.18) and (6.19a) the 
following 







"0 





0" 








"0 





0" 




K = 


^large 











^CKM 




VcKM 























1 














1 



AL (6-24) 



We then obtain 



M: 



33 



(6.25) 



Ko = Vl 



arge 





1 



(6.26) 



where 



1 



M: 



33 



is the 33-matrix element of [M^ ^] , and where we have used the 



CKM properties in the product in (6.24), namely small 

gles, implying Cj w 1. The rotation matrix Viarge is the CKM-analog [10] given 

by 



^larqe — 



C2C3 
— C1S3 — S1S2C3 
S1S3 — C1S2C3 



C2S3 
C1C3 — S1S2S3 

—SiC^ — C1S2S3 



S2 ' 

SIC2 
CIC2 



(6.27) 



where, in order to distinguish it from the CKM matrix, we have used capital let- 
ters. 

Therefore, the neutrino mass matrix, kq is given by 



■ SI 

C2S1S2 



C2S1S2 C1C2S2 

^2*^1 O1O2O1 



O1O2O2 Oi02>Jl ^1^2 



(6.28) 



which is similar to U except for the presence of large angles, and the seesaw con- 
tribution. 
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Needless to say, the parameters of Viarge in (6.27) and (6.28) are what we 
have defined as "input" parameters. 



VII. Mass eigenvalues and mixing angles in the lepton sector with 

texture zeros 

For the charged leptons, as we mentioned earher, the structure of Eq is 



assumed to be the same as that of Dq 











X 





Eq — 


X 


X 


X 







X 


X 



for Ko, as stated in the Introduction, we have two types of structures[18]. 



(a) A-type structure [18] 





X 



X 
X X 
X X 



or kq 



0X0 
XXX 
Q X X 



(7.2) 



matrices give hierarchical neutrino mass values. We will only consider the first 
matrix above since the second matrix gives the same results [18,19]. 
(b) C-type structure [18] 



X 


X 


X 


X 





X 


X 


X 






(7.3) 



gives an inverted hierarchy [18,19]. 

In the following we will not consider the C-typc structure for a simple rea- 
son. Our entire construction of the Yukawa matrices has had at its basis the 
"hierarchical" and " semi- hierarchical" primordial systems which clearly can not 
be compatible with inverted hierarchy. A quick calculation confirms this con- 
clusion. 

We will revisit the C-type structure in section VIII. 

A. Mass eigenvalues for charged leptons 

We simply follow the results from the D-matrix, and write 



En 




e's3 




1 



(7.4) 



the angles s, are the usual "input" CKM angles, and e' is the analog of e 
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in Dq which is the mixing parameter between the "hierarchical" and "semi- 
hierarchical" representations. For the "input" CKM values given in (5.10) and 
e' = .02, we find 



m 

m, 
m 



2. 1 X 10-2 (exp w 6 X 10"^) 
^=9.2x10-4 (exp« 3x10-4) 



(7.5a) 
(7.5b) 



"exp" means experimental values [22]. Once again, there is a good agreement 
with experiments. 

B. Mixing angles and mass eigenvalues for Majorana neutrinos 

The A-type Majorana matrix, kq, has texture zeros as given by [18,19]. 



Ko = 





cqsj 

Cl Cn S 1 



C1C2S2 



clcl 



(7.6) 






CIC2S2 

where Cj and Si are the parameters of the "input" CKM-analog, already defined 
in (6.27) 



y(m) 
large 



C2C3 C2S:i S2 

— C1S3 — S1S2C3 C1C3 — S1S2S3 S1C2 (7-7) 
S1S3 — C1S2C3 —S1C3 — C1S2S3 C1C2 

This matrix is defined in the basis when Eq is diagonal. The experimental values 
of the parameters are [3,4,5,6,7,8] 



5i = (.54 - .83) ; S3 = (.40 - .70) ; ^2 < .16 



(7.8) 



The matrix that will diagonalize kq in (6.28) is 



^'2^3 



8^83 - 



C182C3 



C[C'3- 
— S1C3 



818283 

- C 1 828 3 



S'2 ■ 

8[C'2 
C[C'2 



(7.9) 



However, this matrix is not exactly the CKM-analog [10], V^i°g*e - We need to 
include rotations that diagonalize i?o as well, since the CKM-analog is defined in 
the mass basis of Eq. We note, however, that the rotation angles to diagonalize 
the 2-3 and 1-3 sub matrices of Eq are very small and can be ignored. We will 
only consider the 1-2 submatrix in Eq since, here, for e' of the order .02, needed 
to give the correct charged lepton masses, as discussed earlier, the S3 parameter 
turns out to be w .22 [22]. Thus we define 



V, 



[ouL) 



large 



^2^3 


^'283 


S'2 ' 




C3 


-S3 





r^t Qt Ql Q/ r^l 
— — 0]^02*-^3 


r^l r^l Qf Qf Qf 
^l'-^3 '-'l'-'2'-'3 


8[C'2 




S3 


C3 





Ql Ql r^i Ql r^i 
_ '-'I'-'S ~ <-^l'-'2'-^3 


Ql r^l /^l Ql Ql 

— — LJ1D203 


CiC2_ 










1 



(7.10) 



23 



The effect of tlie seeond bracket above is to change only the angle 9'^, by 
where sinO^E = S3 (~ -22). Thus 



a(out) _ ni a(out) _ ni a(out) _ ni n 
a I — (fl, 02 — (72, (73 — (73 — t/sB 



(7.11) 



First of all a simple examination of the structure of the 23-sub matrix of 
kq in (7.6) shows that the angle of rotation, 6[, that would diagonalize kq is 
precisely ^1 itself, 



(7.12) 



actual value we take for Si will be determined when we calculate the neutrino 
masses. 

For the other two angles we find 



^'2 ~ C1S2 



(7.13) 



where we take C2 = C2 = 1 since ^2 is small. And for S'^,we obtain the following 
equation 



- 5*2 ( ^ ) - 1 = 



(7.14) 



T3 — tan9'^ = 

The diagonalized form for the full mass matrix k, defined in terms of kq in 
(6.25) and (7.6), is found to be the following (again we take C2 = C2 = 1) 



M33 



C*itS'i>S'2 

T' 
-'a 














— T3CiiS'i5'2 
1 



(7.15) 



We now compare the above mass values with the neutrino oscillation data 
which give the following mass-squared differences[3,4,5,6,7,8] 



Am^2 = (2 - 50) X IQ-^el/^ Am^3 = (1.2 - 5) x lO'^eF^ 



(7.16) 



The mixing angle values that are consistent with the experimental values 
given in (7.8) and (7.16) are found to be 



5-2 = .16, Si = .54 
With these parameters and equation (7.14) we obtain 



53 = 0.74 



(7.17) 
(7.18) 
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and from the relation between 0^'^^', ^3 and 63E given in (7.11) we obtain 



4°"*^ = .57 (7.19) 

The mass values for kq are determined from the above parameters and by 
taking 



2 

= 7.1 X IQ-^eV (7.20) 



M33 



where ms^ is the dominant 33-component of the (Dirac) neutrino Yukawa ma- 
trix, N. We then obtain 

mi = 4. 8 X 10-^ey (7.21a) 
m2 = 5. 8 X 10-^ey (7.21b) 

ma = 7.1 X IQ-'^eV (7.21c) 

These values are consisitent with the experimental results (7.8). 

As for comparing the "input" and "output" values of the mixing matrices, 
we have 

g(out) ^ g{in) ^ 54 (7_22a) 

giout) ^ g(in) ^ ^^_22b) 

Wc point out that through (7.9) and (7.10) we first obtained S'g""*-' and then 

simply assumed S'g*"' to have the same value since 5*3 ^= S'g*"'^ was not involved 

in constructing the mass matrix (7.6) for kq and, therefore, there were no con- 
straints on its value. 

The above predictions are consistent with experiments. Moreover, apart 
from the differences in 5*2, which are minor, the "input" and "output" mixing 
parameters are the same. 

We also point out that, in order to get the correct neutrino masses, the 
value of 5*2 needs to be large, as already noted for the case of A-type structure 
in ref.[18] and [19]. As for the possibility of having a maximal coupling in the 
23 sector of ko (i.e. •^'i = = -71) we opted instead to reconcile the "input" 
and "output" values of 5*2 which implied, through relation (7.13), that we have 
as large a Ci as possible, within experimental bounds, which led us to the value 
of 5i given by (7.17). 
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Wc can estimate the mass, M33, for the seesaw neutrino in (7.20) by following 
our assumption that N and E have mass properties that are similar to U and 
D respectively, in which case 

^ « (7.23) 
rriT mh 

and, therefore, putting in the other known mass values, we get m^^ w IQ^GeV'^. 
From (7.20) we then find 

M33 w IQ^^GeV (7.24) 



VIII. Inverted Hierarchy in Neutrinos, revisited 

One might ask if in the eigenvalue equation for the leptons (6.9a,b,c), we 

could have chosen the following "primordial" solutions to generate an inverted 
hierarchy for the neutrinos, leaving the charged leptons in the hierarchical pat- 
tern 

Ai = 1, A2 = 0, A3 = ^{neutrinos) (8.1a) 



Ai =0, A2 = 0, A3 = 1 



or Ai = 0, A2 = 1, A3 = 1 {charged leptons) 

(8.1b) 



First of all, we notice immediately that for the charged lepton RGE (6.1) 
the coupling term 



1 

12' 



NN+E 



will vanish as it involves the product 



"1 





0" 




"0 





0" 







































1 



(8.2) 



(8.3) 



In the absence of the coupling contribution, the charged lepton matrix will have 
two choices, instead of the previous situation when it was, like Dq, a mixture of 
"hierarchical" and " semi- hierarchical" bases. Either it has to be strictly "semi- 
hierarchical" , which will be inconsistent with experiments since it would mean 
two mass eigenvalues very close to each other. Or it has to be "hierarchical" 
which would mean it would be like the up-quark matrix, Uq, which would, 
however, give mass ratios an order of magnitude smaller than observed. Thus 
this model will not work for the charged leptons. We will proceed anyway, with 
the hope that there may be another, totally different, scenario possible for it, 
and take the charged lepton matrix to be diagonal so its parameters do not 
enter into the calculation for neutrinos. 
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The neutrino mass matrix kq will have the following expression (before incor- 
porating texture zeros) 



Ko = Via 



1 





Mn 



(8.5) 



the 11-element of 



with Viarge 8is the "input" matrix as before, and 

Instead of writing a long and complicated matrix that will result from the 
above product we write it after the texture zeros, of the C-variety, are already 
imposed. At the same time, we take account of the fact, evident from kq's 
structure given in (7.3), that with vanishing 22 and 33 diagonal elements, the 
rotation angle for the matrix diagonalizing the 23-sub matrix will be 45*^. Thus, 
the "output" Si will automatically be 



5i 



(out) 



1 

72 



(8.6) 



the above value is, indeed, allowed by the experiments we also take 



1 



(8.7) 



further simplication occurs because of (8.6), namely [19] 

4°"*^ = 

and, here again, we take the same value for the "input" parameter 



(8.8) 



S. 



(in) 







(8.9) 



since it is allowed by the experiments. 

Only ^3 now remains to be determined. With the above values of the "input" 
parameters we obtain the following, simplified expression. 



' CI 


C3S3 
V2 


C3S3 





V2 


C3S3 


02 


I V2 


2 



02 
•^3 



i.lO) 
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Diagonalizing it we obtain ^with T3 = 



rpl 














{-nc^s^ + \si) 

- (I) S: 



(8.11) 



where T3 satisfies 



r<i 1 c2 

'-^S 2*^3 



1 = 



(8.12) 



We choose 53 (= S'-^ = .58, which is within the experimental range (7.8). 
From equation (8.12) we obtain, 53 = .51 which is also the "output" value, since 
we have taken Q-^e = 0, as explained above. Thus we have 



g{out) ^(i") 



V2 



(8.13a) 
(8.13b) 



~i{out) 



.51, 



-,{in) 



.58 



(8.13c) 



The "output" and "input" values are, therefore, identical for two angles and 
very close for the third. 

With the above parameters and with 



Mil 



7.3 X IQ-'^eV 



(8.14) 



before, mfj, 
eigenvalues 



and Mn w 10" GeV) we have the following mass 



mi = 7. 1 X lO-^ey 



(8.15a) 



m2 = 1.2 X 10" W 



(8.15b) 



ms = 8. 3 X 10-^ey 



(8.15c) 



values are consistent with the experimental values given by (7. 16) [3,4,5,6,7,8] 
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IX. Conclusion 

One of our central assumptions was that the scale dependence of a Yukawa 
matrix is dictated entirely by the dominant 33-matrix element which can be 
factored out leaving behind a matrix which is independent of the scale. As 
a conseqence, the renormalization group equations for the Yukawa matrix can 
be expressed as two separate equations: one a differential equation for the 33- 
elcmcnt and another an algebraic equation for the scalc-indcpcndcnt 3x3 matrix. 

It is the properties of the scale-independent matrices that has concerned us 
primarily. After constructing the solutions in terms of the mixing angles and 
incorporating texture zeros, consistent with hierarchical behavior, we made the 
following identifications for U, D, E, and (A-type) k. 





"0 





S2 






' 


es3 0" 




u = 





si 





At, 


D = 




sj + e si 






S2 





1 









Sl 1 





E 








Sl 





Sl 

1 



At 5 






CIC2S2 







CIC2S2 
Ci C2 S 1 



/~< /^2 a /~t2/~f2 
O1O2O1 ^\^2 



M33 



For the C-type neutrinos we have, as discussed in section VII, an especially 

simple mass matrix because of the location of the texture zeros which, automat- 
ically, result in = and 52 = 0, and give 



C3S3 



V2 
C3S3 



C3S3 C3S3 




V2 
Si 



Si 

— # 



^2 



which predicts an inverted hierarchy. 

The above results provide simple expressions for the quark and lepton mass 
matrices in terms of the mixing angles. The manner in which we introduced the 
mixing parameters is self-consistent i.e. what we put in to construct the mass 
matrices (the " input" ) is recovered when we try to obtain the mixing matrices 
(the "output"). 

Another interesting result is that for the massive seesaw neutrinos only one 
mass-scale appears, M33 for normal hierarchy and Mn for inverted hierarchy. 
So the details of the seesaw mass distribution do not play a role, which is an 
important advantage in the model. 

The texture zeros play a crucial role in determining the physical parameters 
in our model. It is well known that zeros in the mass matrices can arise through 
discrete symmetries [1, 11, 23, 24, 25, 26, 27, 28] . The precise group structure we 
need in order to reproduce the above texture zeros needs to be worked out and 
is being attempted currently [29]. 
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To conclude, wc have succeeded in obtaining a simple, transparent, and 
uniform framework to descibe four different pieces of data involving quarks 
and leptons. We assumed the Yukawa matrices to be real and symmetric with 
the dominant eigenvalue as an input. No attempt was made to do a detailed 
numerical analysis to fit the data but the matrices described above are found to 
give a very good descritipn of the mass eigenvalues and mixing matrices. Apart 
from the seesaw particles that are essential to the neutrinos, no new particles 
have been proposed-that is another big advantage to our model. 
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Appendix I 

We notice from (3.9) and (3.10) that, because the determinant and trace are 
invariant, we have 

detM^'^ = detMj^l = (A.l) 



r(i) 



(A.2) 



for simplicity we will take Mq ' to be real and symmetric. Using (A.2) in equa- 
tion (3.4) we obtain 



M, 



(1) 



(A.3) 
which 



we will rewrite as 



l(l+Mr)(l-Mr)MoW=0 (A.4) 

Out of the three factors above, only ^l+Mg^'^has an inverse since its de- 



terminant does not vanish, 



det (l+Mo'') = det (i+M;4) ^ 



(A.5) 



whereas the other two have vanishing determinants. Removing this term in 
(A.4) we obtain, 



(l-M-) M(^) = 



(A.6) 
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Wc note that a general 3x3 real, symmetric matrix with unit trace can be 
written as 



M^'' = {a + b+iy 

And, therefore, 

(l-Mo'') = (a + 6+1)-' 
From (A.6), (A.7) and (A. 8) wc have 



1 + 6 -d -a' 
-d 1 + a -j3 
—a —(3 h + a 



a 


d 


a 


d 


b 





a 


/3 


1 



-a 



1 + 6 -d 

~d 1 + a -/? 
—a —(3 6 + a 



- + a{h+l) 
-ad — aj3 + d (o + 1) 
-aa — djS + a{a + b) 



-d^ -13^ + b 
-bp -da + P{a + b) 



a 


d 


a 


d 


b 





a 




1 


HI) 




—a 


-1) 




-/3 



(A.7) 



(A.8) 



a + b-a^-l3^ 



(A.9) 

= 
(A.10) 



It is then easy to show that 



matrix Mq'^ is then given by 



a = a 
b = 0^ 

d = afi 



M, 



(1) 





a/3 


a 


aft 






a 


f3 


1 



If the 33-term above dominates the matrix so that 



a < 1 and /? < 1 
then we have a classic hierarchy pattern in Mq^\ 



(A.lla) 
(A.llb) 

(A.Uc) 



(A.12) 



(A.13) 
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Appendix II 

Equation (4.22) given by 



dx 2 



12 



(B.l) 



can be solved analytically by first substituting expression (2.3) for and solving 

(B.2) 



dy ^ .2 1 .2 I. . .2 \-i 



which gives, with C = constant, 



If we then take 



and normalize 



2 ^\~T2 



Xb{x)=f3{x){l + Xlx) 



Xb{0) = /3(0) = Xob 



(B.3) 



(B.4) 



(B.5) 



then the equation for f3{x) is obtained from (B.l) whose solution is then given 
by the following complicated relation 



1 6 



5X1, 



{1 + Xl,x) ^-1 



1 



A, 



(B.6) 



05 



Substituting this value of P(x) in equation (B.4) will then give us the com- 
plete analytic expression for A(, 
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